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Abstract 


Space radiation transport codes require accurate models for hadron production 
in intermediate energy nucleus - nucleus collisions. Codes require cross sec- 
tions to be written in terms of lab frame variables and it is important to be able 
to verify models against experimental data in the lab frame. Several models 
are compared to lab frame data. It is found that models based on algebraic 
parameterizations are unable to describe intermediate energy differential cross 
section data. However, simple thermal model parameterizations, when appro- 
priately transformed from the center of momentum to the lab frame, are able 
to account for the data. 


1 Introduction 

Many radiation transport codes are used for transport of nucleons and other hadrons, 
and are suitable for simulating such things as particle physics experiments and cosmic ray 
proton propagation through interstellar space and the Earth atmosphere. These codes 
are also used in proton and hadron therapy applications. Some codes can transport heavy 
ions and can be used in space radiation applications as well as heavy ion beam therapy. 
An example is the code HZETRN [1, 2, 3] which is a heavy ion transport code widely 
used in space applications. Heavy ion transport codes need many atomic, nuclear, and 
particle physics cross sections as input. The atomic cross sections are needed to describe 
slowing down of the incoming projectile. A large range of nuclear interactions cross 
sections including fragmentation, breakup, electromagnetic dissociation [4], and even pair 
production [5] are needed as input. 

A special difficulty of space radiation transport codes is the intermediate energy region 
in which they must be accurate. Low energy cross sections can be dealt with using 
non - relativistic formalisms and very high energy cross sections can also be dealt with 
using ultrarelativistic approximations. However, the peak of the cosmic ray spectrum [6] 
occurs in the intermediate energy range of several GeV. From a theoretical view, this is 
a difficult area because non - relativisitic theory is not applicable and ultrarelativistic 
approximations do not work. Eyser and Machleidt [7] have emphasized this with regard 
to nucleon - nucleon scattering. Heavy ion transport codes also need to include hadron 
production from nuclear collisions, which is the subject of the present paper. 

Another difficulty faced by many transport codes, including heavy ion space radiation 
codes, is that cross sections must be calculated in the lab frame because the Boltzmann 
transport equation is usually set up in such a way that projectiles are transported through 
a stationary target. Particle physics cross sections, however, are much easier to calculate 
in the nucleon - nucleon center of momentum (cm) frame. Although the techniques 
for Lorentz transforming to the lab frame are well known, one must deal with extra 
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complications such as double valued functions of angles [8] and the fact that cross sections 
in the lab frame are much more complicated than their relatively simple counterparts in 
the cm frame. In summary, some of the difficulties faced by heavy ion transport codes for 
space radiation applications are: 

1) They must include all relevant atomic, nuclear and particle physics effects. 

2) They must be accurate in the theoretically difficult intermediate energy range. 

3) All cross sections must be transformed to the lab frame. 

The aim of the present paper is to discuss parameterizations of pion and kaon produc- 
tion in heavy ion collisions at intermediate energy in the lab frame. These parameteriza- 
tions will be suitable for use in heavy ion transport codes for space radiation applications. 


2 Differential cross sections 

The Lorentz - invariant differential cross sections of Badhwar [9] and the thermal model 
have very simple forms and they are both expressed in the nucleon - nucleon center of 
momentum frame. Even though they are Lorentz - invariant, we shall always want to plot 
them as functions of non - invariant variables and indeed they are written in terms of non 
- invariant variables. The Nagamiya data [10] is given in the lab frame and so we must 
transform the simple looking Badhwar and thermal model cross sections into complicated 
lab functions. Furthermore, transport codes require all cross sections in the lab frame. 
The techniques for performing the Lorentz transformations are described below. 

2.1 Review of kinematics 

The appropriate relativistic kinematics is reviewed in this section, so that one can under- 
stand the quantities used in the algebraic parameterizations, which are introduced later. 
Consider the inclusive reaction 


a T b — > c T X , (1) 

where c is the produced particle of interest and X is anything. Throughout this paper we 
assume that all variables, such as all momenta, are evaluated in the center of momentum 
(cm) frame, unless otherwise indicated. Lab frame variables will be given a subscript. For 
example, the variable x evaluated in the cm frame is written as x but evaluated in the lab 
frame, it is written as x i a b- The momentum of particle c is denoted as p, and supposing 
that it scatters at angle 6 to the beam direction, then the longitudinal and transverse 
components of momentum are 


= pcosO , 
= psind . 


Pz 

Pt 


2 


( 2 ) 

( 3 ) 



Note that 


2 2,2 
P = Pz+Pt 1 

tan 9 = Pt/Pz ■ 


(4) 

(5) 


Feynman used a scaled variable instead of p z itself [11, 12, 13, 14], The Feynman scaling 
variable is [13, 15, 16, 17, 18, 19, 20] 

x F = , (6) 

Pz max 

where p z is the longitudinal momentum of the produced meson in the cm frame, and 
Pz max is the maximum transferable momentum given by [16, 17, 20] 


Pz max 


l\(s,m c ,m x ) 
4s 


(7) 


with 


A(s, m-i, rrij ) = (s — rri 2 — m 2 ) 2 — 4 mfm 2 . (8) 

Here, the Mandelstam variable s is the square root of the center of momentum total 
energy, and m is the mass of a particle. Note that 

Pz max = Pmax ) (9) 

because the component of any vector can be no larger than the magnitude of the vector 
itself. Nagamiya and Gyulassy [17] point out that if c is a boson with zero baryon number, 
then 


rn x = m A + m B , 


( 10 ) 


in agreement with the p z max formulas of Nagamiya and Gyulassy [17] and Cassing [20]. 
The Feynman scaling variable approaches the limiting value [18] 


Xf 



as 


oo 


( 11 ) 


Also, it is obviously bounded in the following manner [13] 


— 1 < Xp < 1 • 


( 12 ) 


Sets of variables that are often used are either (p, 6) or ( p z ,pr )• Writing 


Pz 


Xp 


' X(s, m c , m x ) 
As 


(13) 
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shows that another useful variable set is (. x F ,p F ), which is used by Alt et al. [21, 22] when 
presenting their data. These variables are also used throughout the present work. 
Rapidity is defined as 


so that 



f E + p z \ 
\E — p z ) 
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E = niT cosh y , 

p z = TYIt skill y , 

where the transverse mass is defined through 

2 2 i 2 771 2 2 

m T = m + p T = E — p z , 


(14) 

(15) 

(16) 

(17) 


with m is the mass of the produced particle c. This gives yet another useful variable 
set (; u,Pt )• In the work below, it will be necessary to write the rapidity in terms of the 
Feynman scaling variable as 



yJx 2 F + m%/iP zmax + x F ' 
lx 2 F + m 2 T /p 2 zm ^-x F/ 


(18) 


2.2 Lorentz transformations 

The relativistic gamma factor for the projectile nucleus in the lab frame is given by 

7p = T p /m + 1 , (19) 

where T p is the kinetic energy of the projectile nucleus in units of energy /nucleon and m 
is the mass of the nucleon. Using 

gives the speed (3 P of the nucleus in the lab frame. This is also the gamma factor for a 
particular nucleon in the lab frame and f3 p is the speed of a particular nucleon in the lab 
frame. 

Denote the nucleon - nucleon center of momentum frame as NNcrn and the nucleus - 
nucleus center of momentum frame as AAcm. We want the speed of the NNcm in the lab 
frame, which we denote simply as (3. The relativistic gamma factor for the NNcm in the 
lab frame is 

7 = ’ 



( 20 ) 
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The relation between the two gamma factors is given by 

Tv = 27 2 - 1 • (22) 

This relates the speed of the NNcm in the lab frame to the speed of a projectile nucleon in 
the lab frame. The gamma factor that appears in the Lorentz transformations below is 

7 = \J\( 7„ + l) ■ (23) 

This 7 enables us to transform quantities from the NNcm frame to the lab frame. One 
can now see how this is obtained from the nucleus kinetic energy T p in equation (19). 
The cm frame moves at speed f3 relative to the lab. The Lorentz transformations are 

E \ = / 7 -7/5 W ^lab \ 

Pz ) V 7 / \ Pz lab j 

with the transverse momentum p t being the same in both frames. The parameteriza- 
tions below are often written in terms of the cm variables ( Pz,Pt )• The explicit Lorentz 


transformation to (p z i a b,PT lab) variables is given by 

Pz = ~TP\JpI lab + Pt lab + + 7 Pz lab , (25) 

PT = PT lab • (26) 

Instead of (p z i a b, Pt lab) variables, we will also need to use (pi a b, $iab) variables. The explicit 
transformation to the momentum and angle lab variables is 

Pz = “7/5 \/pL + m2 + 7Piab cos 0i ab , (27) 

Pt = piab sin 6 > lab • (28) 


Simple substitution enables one to write the rapidity and Feynman scaling variable in 
terms of lab variables. 

2.3 Algebraic parameterizations 

Badhwar and Stephens [9, 23] have provided parameterizations of inclusive pion and kaon 
production from proton - proton collisions. These parameterizations are for high energy 
and are given in the nucleon - nucleon center of momentum frame (NNcm), and they are 
given for both charged and neutral pions and for charged kaons. These parameterizations 
have been transformed to the lab frame using the techniques above. The results for 
charged and neutral pions are plotted in figure 1 . As expected, the cross sections for 7 r + 
are the largest, tt~ are the smallest and n° are in between. The charged kaon cross sections 
are plotted in figure 2. It can be seen that the parameterizations are much smaller than 
the pion cross sections and again the K + cross sections are bigger than the K~ cross 
sections. Other parameterizations have also been developed by Alper [24], Mokhov [25], 
Ellis [26] and Carey [27]. Details of the various parameterizations are now given. 
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2.3.1 Badhwar parameterization 

The Badhwar parameterization [9] gives the Lorentz - invariant differential cross section 
for charged pions as 


d 3 p 



and neutral pions as 


A 

(1 + 4 m? p /s) r 


(1 — :r) 9 exp[- 


-Bpr -I 
+ 4 m 2 p /s 


d 3 p 



Af(E p )(l - x) q ex p 


- Bp T 
1 + 4 rn^/s 


1 


and charged kaons as 


d 3 p 


(K±) 


A(1 — x) c exp (—Bp T ) , 


(29) 


(30) 


(31) 


where m p is the proton mass, y/s is the total energy in the center of momentum (cm) 
frame, and pr is the transverse momentum of the produced meson in the cm frame. The 
other terms are given by 


x 


x, 


4. ■ 

+ - (z* 

s 


+ m 2 ) 


I 1/2 


(32) 


where it is assumed that the variables appearing in Xp are in the cm frame. The mass m 
is the mass of the produced particle (pion or kaon). Badhwar writes xt = xf- Also, 


C\ + C-2Pt + Cs pp 
y^l + Amy s 


(33) 


The function f(E p ) for neutral pions is given by 

f(Ep) = (1 + 23iT- 2 - 6 )(l - 4 ml/sY , (34) 

with the constants listed in Table 1. Badhwar points out that for large values of E p , 
equation (30) takes the asymptotic form 

d 3 n 

EEp-(7T 0 ) = Aexpi-BprMl-x)^- 0 ^ 0 ^ , (35) 

d 3 p 

which is consistent with the Feynman scaling hypothesis [13]. The Badhwar variables are 
{xpi Pt) i which are also used in the Alt et al. [21, 22] data, so that no variable conversion 
is necessary. 
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Table 1: Constants for the Badhwar parameterization. 


Particle 

A 

B 

r 

C 

Ci a 

c 3 

7T + 

153 

5.55 

1 


5.3667 - 3.5 

0.8334 

7 

127 

5.3 

3 


7.0334 - 4.5 

1.667 

7T° 

140 

5.43 

2 


6.1 3.3 

0.6 

K+ 

8.85 

4.05 


2.5 



K~ 

9.3 

3.8 


8.3 




2.3.2 Alper parameterization 

The Alper [24] parameterization for charged pions and kaons, protons and antiprotons is 

7 2 (1 - p T /Pbeam) m 

E lfp = A exp ( r Bpr) exp {-Dy ) + A 2 ^ + ^ , (36) 

with the constants listed in Table 2. The Alper variables are ( y,PT )• To change to the 
variables ( xp^Pr ), we convert the rapidity in equation (36) to xp using equation (18). 


Table 2: Constants for the Alper parameterization. 


Particle 

A\ 

B 

D 

A 2 

M 

m 

n 

7T + 

210 

7.58 

0.20 

10.7 

1.03 

10.9 

4.0 

7T~ 

205 

7.44 

0.21 

12.8 

1.08 

13.1 

4.0 

K+ 

14.3 

6.78 

1.5 

8.0 

1.29 

12.1 

4.0 

K~ 

13.4 

6.51 

1.8 

9.8 

1.39 

17.4 

4.0 

P 

5.3 

3.8 

- 0.2 

16 

1.2 

0 

7.5 

P 

1.89 

4.1 

2.3 

25 

1.41 

25 

4.5 


2.3.3 Ellis parameterization 

The Ellis [26] parameterization for charged pions, neutral pions, charged kaons, protons 
and antiprotons is 


d 3 p 


A(p 2 T + M 2 y N/2 ( 1-xt) F , 


(37) 
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where A is an overall normalization fitted to be A = 13 in reference [28] and Xp = 
Pt / Pmax ~ 2 'pr/y/s. The same value of A is used in the present work. The other con- 
stants are listed in Table 3. The Ellis parameterization is independent of the emission 
angle 6, and so does not carry any dependence on p z , xp, y etc. 


Table 3: Constants for the Ellis parameterization. 


Particle 

N 

M 2 

F 

7 r + 

7.70 

0.74 

11.0 

7T~ 

7.78 

0.79 

11.9 

7T° 

10.8 

2.3 

7.1 

K+ 

8.72 

1.69 

9.0 

K~ 

8.76 

1.77 

12.2 

P 

10.38 

1.82 

7.3 

P 

9.1 

1.17 

14.0 


2.3.4 Mokhov parameterization 

The Mokhov [25] parameterization is 

E A^ = A ( 1 P ) ex P V i(Pt)V2(pt) , (38) 

d A p \ Pmax J \ Cy/Sj 

where 

Vi (pt) = (1 — D) exp(— Eptp) + D exp (— Fp^) for p T < 0.933 GeV, 

0.2625 

Vi (pr) = ^2 + Q 87)4 for p T > 0.933 GeV, (39) 

and 

V 2 (p T ) = 0.7363 exp(0.875pT) for pp < 0.35 GeV, 

V 2 (p T ) = 1 for p T > 0.35 GeV, (40) 

with the constants listed in Table 4. Using p = \Jp 2 z + Pp, gives the Mokhov variables 
( PziPt ) which are transformed to ( xp,Pr ) using equation (13). 



Tabic 4: Constants for the Mokhov parameterization. 


Particle 

A 

B 

C 

D 

E 

F 

7 T + 

60.1 

1.9 

0.18 

0.3 

12 

2.7 

7T~ 

51.2 

2.6 

0.17 

0.3 

12 

2.7 


2.3.5 Carey parameterization 

The Carey [27] parameterization, for negative pions, negative kaons, and antiprotons is 

rjZfj 

E lfp = hN (PT + G)~ 4 - 5 (l -x R ) J 1 (41) 

where N is an overall normalization fitted to be IV = 13 in reference [28] and Xr = 
p/Pmax ~ 2 pl\[s. The same value of N is used in the present work. The constants are 
listed in Table 5. The Carey variables are ( P z ,Pt )• To change to the variables (x f ,Pt), 
we use x R = \jx 2 F +p 2 T /p 2 m ^. 


Table 5: Constants for the Carey parameterization. 


Particle 

N 

h 

G 

J 

7 

13 

1.0 

0.86 

4 

K 

13 

0.36 

1.22 

5 

P 

13 

0.26 

1.04 

7 


2.4 Thermal model 

Several authors [10, 17, 29, 30] have shown that a simple thermal model parameterization 
is able to describe both pion and kaon data in intermediate energy nucleus - nucleus 
collisions. In the nucleon - nucleon cm frame, this parameterization is given by 

E^- = NeM~T/E 0 ) , (42) 

where T is the meson kinetic energy in the nucleon - nucleon cm frame. (It is not the 
nucleus - nucleus cm frame.) The parameters N and Eq are constants fitted to the data. 


9 



These parameters are different for protons, pions and kaons. Some effective nuclear tem- 
peratures Eq are listed in Table 6, which have been taken from references [10, 17, 30]. 
We use the wording “effective” nuclear temperature instead of just nuclear temperature 
because there is not a single well defined temperature. It is different for different particles, 
such as pions and kaons. Equation (42) is characteristic of a thermal Boltzmann distri- 
bution, and describes the “boiling” off of hadrons after nuclear equilibration has been 
reached. Clearly, this model is vastly different from a nucleon - nucleon collision model 
assumed above. 


Table 6: Effective nuclear Boltzmann temperatures Eq. 


Initial State 

Projectile 

Produced 

Eq 

References 


Energy 
(A GeV) 

Particle 

(MeV) 


Ar + Pb 

0.8 

7T~ 

72 

figure 19 [10] 

Ar + Pb 

0.8 

P 

87 

figure 19 [10] 

Ar + KC1 

0.8 

7T~ 

66 

figure 19 [10] 

Ar + KC1 

0.8 

P 

79 

figure 19 [10], figure 4 [30] 

C + C 

0.8 

7T~ 

60 

figure 19 [10] 

C + C 

2.1 

7T~ 

123 

figure 27 [17] 

C + C 

0.8 

P 

68 

figure 19 [10], figure 4 [30] 

Ne + NaF 

0.18 

7T~ 

35 

figure 27, 28 [17] 

Ne + NaF 

0.4 

7T~ 

39 

figure 28 [17] 

Ne + NaF 

0.8 

7T~ 

62 

figure 19 [10], figure 28 [17] 

Ne + NaF 

2.1 

7T~ 

102 

figure 28 [17] 

Ne + NaF 

0.8 

P 

75 

figure 19 [10], figure 4 [30] 

Ne + NaF 

2.1 

K+ 

142 

figure 35 [17], figure 5 [30] 
See also [29] 


3 Comparison to experiment 

Experimental data for it ~ production in the Ar + KC1 nuclear collisions at 800 A MeV have 
been measured by Nagamiya et al. [10]. They measured Lorentz - invariant differential 
cross sections as a function of transverse momentum and lab angle. 

The nuclear radius is given by the well known relation Rth 1.2 fin A 1 / 3 , and therefore a 
total cross section for a nucleus scales like A 2 / 3 multiplied by a nucleon cross section. The 
nucleus - nucleus geometrical cross section [10, 17, 31] therefore scales as (Ap 3 + A^ 3 ). 
However, pion production is more dependent on the number of binary nucleon - nucleon 
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collisions, which has a different scaling function [10]. The Badhwar parameterization 
is scaled up to nucleus - nucleus collisions by multiplying the nucleon - nucleon cross 
section by {ApAp ) 2 ^ , where Ap and Ap is the projectile and target nucleon numbers. As 
discussed above, this is then transformed to the lab frame. The results are compared to the 
Nagamiya data in figures 3 and 4. It can be seen that the scaled Badhwar parameterization 
fails to describe the data. Even if the scaling factor is changed to an arbitrary constant, 
the agreement with the spectral shapes is still very poor. Extensive work has also been 
done on the parameterizations of Alper [24], Mokhov [25], Ellis [26] and Carey [27]. These 
parameterizations have been treated in the same fashion as the Badhwar parameterization 
and compared to data. All of them also fail to describe the data, with comparisons of 
similar quality to the Badhwar comparison. 

The failure of algebraic parameterizations is not unexpected. They were developed to 
describe high energy nucleon - nucleon collisions in the ultrarelativistic limit, whereas the 
Nagamiya data describe intermediate energy nucleus - nucleus collisions. 

The pion thermal model parameterization, transformed to the lab frame, is shown in 
figures 5 and 6. The parameter N is fitted to the data with 


N 


1.7 x 10 4 


mb 

GeV 2 sr 


0.017 


mb 

MeV 2 sr 


The effective nuclear temperature is taken from Table 6 which is Eq = 66 MeV. Com- 
parison to experiment is very good. It is clear that the rather complicated experimental 
distributions as a function of lab angle reveal themselves to be nothing more than a simple 
NNcm thermal distribution appropriately transformed to the lab frame. Finally, the kaon 
thermal model parameterization, transformed to the lab frame, is shown in figure 7. The 
reason this is plotted is for comparison to the pion spectra. 


4 Extension to other nuclei, energy and pion species 


Nagamiya et al. [10] have discussed pion production from arbitrary nuclei at arbitrary 
energy. This is needed for transport codes such as HZETRN. The comparisons presented 
so far have been for the Ar + KC1 system with projectile mass number of Ap = 40 and 
target mass number Ap = 37 where the latter is the average of Ak = 39 and Aci = 35. 
Nagamiya et al. [10] show that a simplified dependence is 

d 3 n 

E d?p K ■ (43) 


Equation (42) is true for A P = 40 and Ap = 37. Equation (43) means that equation (42) 
can be generalized to 


d 3 p 


/ ApAp \ 5/3 

V40 x 37/ 


N exp(—T/E 0 ) 
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for arbitrary values of Ap and Ap. In figure 25 of reference [10], the mass dependence of 
the effective nuclear temperature Eq is shown. This temperature can be parameterized as 


E 0 = 70/r lab /1000 , (45) 

with Ti ab and E 0 both in units of MeV. Thus, the final generalized cross section is 

d 3 a /A P A T \ 5/3 ( T \ 

d 3 p V40 x 37 ) ^ 70^T lab /1000 ) 

with all energies in MeV, resulting in the cross section E^ with units of M ^ 3sr . This 
generalization is expected to be accurate only for energies Ti ab < 5 AGeV and should not 
be used for higher energy. 

Finally, the cross sections developed in the present work have been applied to n~ 
production only. From figure 30 of Nagamiya et al. [10], one can see that to a rough 
approximation , the ratio of positive to negative pions is roughly unity. Employing the 
well known relation a(7r°) « |[<t(7t + ) + a(7r - )] shows that the cross section for all pion 
species is the same if the cross sections for positive and negative pions is the same. Thus, 
in the present work, the approximation that the cross section is the same for all pion 
species is used. This approximation should be improved in future work. 

5 Conclusions 

Space radiation transport codes require accurate models describing hadron production in 
intermediate energy nucleus - nucleus collisions. These models must be written in terms of 
lab frame variables. Several models have been compared to the intermediate energy data 
of Nagamiya [10], which has been presented in terms of lab frame variables. The Erst set 
of models used high energy algebraic parameterizations multiplied by (ApAp) 2 / 3 . (Scaling 
by (. ApAp ) 5//3 also fails for these models.) Comparison between theory and experiment 
was poor. Such models should not be used in 3 - dimensional transport codes. A simple 
thermal model parameterization was found to satisfactorily describe the data. Because 
of this success, it is recommended that such thermal model parameterizations should be 
utilized in 3 - dimensional space radiation transport codes. 
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Figure 1: Pion inclusive Lorentz invariant differential cross sections as a function of lab 
momentum and angle for proton - proton collisions with an incident proton kinetic energy 
°f Fiab = 5 GeV. The solid lines are 7r°, the long dashed lines are n~ and the short dashed 
lines are 7r + . The lab angles from top to bottom are 15°, 20°, 30°, 40°, 60°, 90°, 110° and 
145°. The curves have been generated using the Badhwar and Stephens parameterizations 
[9, 23] and by transforming to the lab frame. 
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Figure 2: Kaon inclusive Lorentz invariant differential cross sections as a function of 
lab momentum and angle for proton - proton collisions with an incident proton kinetic 
energy of Tj a b = 5 GeV. The long dashed lines are K ~ and the short dashed lines are K + . 
The lab angles from top to bottom are 15°, 20°, 30°, 40°, 60°, 90°, 110° and 145°. The 
curves have been generated using the Badhwar and Stephens parameterizations [9] and 
by transforming to the lab frame. 
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Figure 5: Thermal spectrum successfully describes data [10] at 800 A MeV. The overall 
normalization from equation (44) was fitted as N = 1.7 x 10 4 mb/GeV 2 sr. The effective 
nuclear temperature, taken from reference [10], is E 0 = 66 MeV. Lab angles are indicated. 
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Ar + KC1 -> K + + X 



Figure 7: K + spectra at 2.1 A GeV, with the effective nuclear temperature, taken from 
references [17, 30], is Eq = 142 MeV. Overall normalization is arbitrary. Lab angles from 
top to bottom are 15°, 20°, 30°, 40°, 60°, 90°, 110°, 145°. 
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